Coexistence of ferromagnetism and singlet superconductivity via kinetic exchange 
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We propose a novel mechanism for the coexistence of metallic ferromagnetism and singlet super- 
conductivity assuming that the magnetic instability is due to kinetic exchange. Within this scenario, 
the unpaired electrons which contribute to the magnetization have a positive feedback on the gain 
of the kinetic energy in the coexisting phase by undressing the effective mass of the carriers involved 
into the pairing. The evolution of the magnetization and pairing amplitude, and the phase diagram 
are first analyzed for a generic kinetic exchange model and then are determined within a specific 
case with spin dependent bond-charge occupation. 



The problem of the interplay between ferromagnetic 
(FM) and superconducting (SC) long range order has 
been recently attracting new interest due to the dis- 
covery_of superconductivity in ferromagnetic metals, 
UGe2;i], ZrZnzJi, URhGe'ill, and in rutheno-cuprate 
RuSr2RECu208 compounds, with RE=Eu or Gd Pj. 

The investigation of ferromagnetic superconductors 
started by analyzing the case of two interacting sub- 
systems: one formed by localized spins or aligned mag- 
netic impurities which produces the ferromagnetic back- 
ground, the other composed by itinerant electrons which 
gives rise to the superconductivity. Within this frame- 
work, early works [a, IE focused on singlet supercon- 
ductivity in presence of a spin-exchange field, showing 
that it can exist only below a critical value of the mag- 
netic coupling. Hence, with the purpose to increase the 
threshold of the critical spin-exchange, it was suggested a 
finite momentum pairing state coexisting with the ferro- 
magnetic order IE 0|- In this pairing configuration, com- 
monly indicated as the Fulde-Ferrell-Larkin-Ovchinnikov 
state, the subtle balance between the condensate energy 
of the Cooper pairs with a finite center-of-mass momen- 
tum and the Zeeman energy, related to the magnetic mo- 
ments of depaired itinerant electrons, was able to account 
for a superconducting-ferromagnetic (SF) phase in pres- 
ence of spin-exchange higher than the zero-momentum 
pairing state. 

Nevertheless, in the above mentioned materials, a new 
phenomenology seems to arise if compared to the conven- 
tional case of a metal with magnetic impurities. Indeed, 
it has been suggested that (i) ferromagnetism and su- 
perconductivity are cooperative phenomena, (ii) the FM 
state is due to itinerant electrons, (iii) the same electrons 
participate in both the FM and SC order, and (iv) for the 
case of systems with two types of carrier responsible for 
the SC and FM phases separately, there are interesting 
cooperative effects due to competing charge- and spin- 
exchange coupling. Theoretical studies along those direc- 
tions have been recently performed by considering the oc- 
currence of ferromagnetism and superconductivity in ei- 



ther singlet [mil El El and/or triplet El El El chan- 
nel of pairing. In particular, as far as the SC state with 
singlet s-wave pairing is concerned, it has been shown0 
that its coexistence with weak itinerant ferromagnetism 
can be obtained within a single band model, where the 
ferromagnetic order is driven by the same electrons that 
participate in the formation of Cooper pairs. One cru- 
cial aspect of such analysis concerns the stability of the 
SF state. The loss of condensation energy due to the 
depaired electrons which produce a non zero total mag- 
netization, is not compensated by the energy gain of the 
magnetic exchange, thus the SF state turns out always to 
be energetically unfavorable against the non-magnetic SC 
one even if s- and d-wave symmetry or finite-momentum 
pairing state are considered [l^. 

In this letter, we propose a novel mechanism for 
the coexistence of superconductivity and ferromagnetism 
within a single band model. The new and crucial ingre- 
dient is that the metallic ferromagnetism is not due, as 
in the previous studies, to a rigid shift in the positions of 
the majority and minority spin bands (i.e. Stoner model), 
but it is a consequence of a change in the relative band- 
width of electrons with up and down spin polarization. In 
this circumstance, the gain of energy comes from the un- 
dressing of the mass for the majority spins which induces 
a bandwidth enlargement and in turn lowers the kinetic 
energy El- When the pairing interaction is switched on, 
the interplay between the gain in kinetic energy, due to 
the relative change of the majority and minority spin 
bands, and the condensation energy which would tend to 
pair all the electrons becomes crucial for the stability of 
the SF state. It turns out that such a phase is the most 
favorable only if a suitable tuning of the ratio between 
the density of depaired and paired states is reached, so to 
optimize the balance in the kinetic and potential energy 
gain. We will analyze this possibility in two steps: I) the 
case in which the variation between the mass of up and 
down spin electrons is arbitrary modified without refer- 
ring to any microscopic mechanism is firstly considered; 
II) a specific tight-binding model in which off-diagonal 
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Coulomb interactions are responsible for an asymmetric 
bond-charge distribution for each spin channel, is then 
investigated. For both cases, the attractive interaction 
will be assumed of BCS-type. 

Model I. Let us start from a model Hamiltonian which 
contains a local attractive potential and itinerant elec- 
trons with a spin dependent mass: 

Hi = ^ (2w<,)(c,Lcj<,-|-i7.c.)-g^4^4^C4Ci| 



netization have the following form 



E4 



(1) 



where cl^{cic) creates (destroys) an electron with spin a 
at the site i, Wa is a positive term that controls the renor- 
malization of the mass for electrons with spin a (the fac- 
tor 2 being introduced for convenience) , /i is the chemical 
potential, t is the hopping amplitude which defines the 
bare bandwidth, and g is the pairing coupling being ef- 
fective only in a shell of amplitude 2wc around the Fermi 
surface as in the usual BCS theory. It is worth pointing 
out that in this case the mass dressing and undressing 
can be generated both via the coupling to a background 
of spin, as in double-exchange systems, or to dynamical 
processes intrinsically generated by electron correlations. 

The introduction of a pairing amplitude, after the 
mean-field decoupling of the attractive term, brings to 
the following diagonal expression for the Hamiltonian: 

Hj.MF = E {^taW + E^kPlPk) + ^0 (2) 
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(3) 



where = — t exp {ik ■ 6) is the bare dispersion with 
5 being a vector connecting a site to its nearest neigh- 
bors, and A — J2k si'^ki'^-ki) pairing amplitude, 
respectively. The field operators a^, correspond to 
fermionic excitations with quasiparticle dispersion 



a,/3 



(4) 



b — + wi and a — w-^ — wi being the average and half 
the difference of the spin mass renormalization. It is im- 
mediately apparent that the kinetic exchange amplitude 
is proportional to a while the value of b determines the 
order of magnitude of the average kinetic energy. Still, 
the sign of the total magnetization (M) follows that of 
a, while 6 is always positive. Thus, without any loss of 
generality one can focus only on the case with a > 0, the 
other one being just derived by inverting the direction of 
M. 

The equations for the pairing amplitude and the mag- 
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(5) 
(6) 



where n'^'^ are the Fermi distributions in the momentum 
space of the fermionic fields having E'^'^ as energy dis- 
persion relations, respectively, and /is is the Bohr mag- 
neton. Let us now discuss under what conditions the 
system can accommodate a coexistence of SC and FM 
order. The analysis will be restricted to the zero tem- 
perature limit and furthermore we will assume that the 
density of electrons n is self-consistently fixed by a posi- 
tive value of the chemical potential /x. The results can be 
symmetrically extended within the same procedure to the 
negative /i case. Hence, considering that &, a, A,/j, > 
and 6 > a, it is possible to show that = as Ej^ is 
always positive, while — 1 when Afc^_ < Cfc < ^k,+ . 



where Xk.± 



The Xk,± define a 



range on the positive side Cfe, which is different from 
zero only when a and b are linked to the amplitude of 



A and /z in a way to fulfill the relation f > 



If 



the previous inequality holds, single particle states form 
within the gap and contribute to give a finite magne- 
tization. Then, depending on their density the system 
can allow for a state with coexisting SC and FM order. 
In the Hmit of /i ^ (n 1) or a (M ~ 0), 
there are no real solutions for Xk,± so that one ends up 
with the usual BCS-like state with zero magnetization. 
In the other cases, the self-consistent equations for the 
gap amplitude and the number of electrons have to be 
solved to determine the conditions for the existence of 
SF state. For such purpose, it has been used a model 
density of states given by N{e) = l/2?x; (with w — 2t) 
if —w < e < w. We notice that, contrary to the case 
of the Stoner model, a non-trivial ferromagnetic solution 
can be got without the need of a curvature in the den- 
sity of states close to the Fermi level|0. The analysis of 
the solutions is performed by fixing the value of b and n 
and studying how the pairing amplitude, the energy and 
the magnetization are modified by changes of a. Due to 
the formation of unpaired electrons within the gap, the 
value of the effective average mass in the SF state (&s/) 
is modified with respect to the case of zero magnetiza- 
tion (6sc). Hence, there might occur two distinct cases 
depending on the microscopic mechanism which controls 
the relative change of the majority and minority spin 
bandwidth: f ) undressing of the average effective mass, 
i.e. 6b > 0; 2) dressing of the average effective mass, i.e. 
Sb <0, 5b= {bsf — bsc)/bsc being the relative percentage 
variation. 

In Figgis shown the comparison between the energy 
of the SF and the SC state, together with the behav- 
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ior of the magnetization and the pairing amphtude as a 
function of the relative shift between the mass of up and 
down spin electrons. For simplicity, we have fixed b^c = 1 
and we have studied how the difference of energy is mod- 
ified in the cases (1) and (2). Changes in bsc produce 
only quantitative but not qualitative differences in the 
results. As shown in the top panel of Fig^ only when 
db is positive and for a suitable density of the depaired 
states, controlled by the amplitude of a, the system can 
be stabilized in a SF phase. When the average mass in- 
creases, though the effective amplitude pairing grows and 
in turn lowers the condensation energy, the loss in the ki- 
netic energy of the electrons that participate in the pair- 
ing cannot be counterbalanced to get a stable SF state. 
It is worth pointing out that the transition from the SF 
state to the SC one is of first order for the magnetization 
but it is continuous for the pairing amplitude by mov- 
ing from small to large values of a, while it is reversed 
in the opposite a direction. Nevertheless, depending on 
the amplitude of the undressing Sb, the transition from 
the SF state can be first order type in both the order 
parameters. Finally, due to the peculiar link between 
A and fi, we notice that the change in the density of 
electrons modifies the region where the SF solution ex- 
ists. In particular, the interval of solutions shrinks if one 
moves towards the half-filling case, whereas a large band 
undressing is required to stabilize the SF phase. 

Model II. Now we consider an explicit tight-binding 
model where the bandwidth change depends on the bond- 
charge occupation in each spin channel [lij. The model 
Hamiltonian has the following form: 

Hii = -t {4aC3<r + H.cj +UY^ rii^nn 

<ij>(T i 

+ V n,nj +J YZ i4aC]a-Cra'Cja + H .c) 

<ij> <ij ><7<7' 

<ij> i 

where U and V represent the on-site and nearest- 
neighbor Coulomb repulsion, respectively, and the pa- 
rameters J and J' describe nearest-neighbor exchange 
and pair hopping processes, while g is the BCS pairing. 
By applying the Hartree-Fock decoupling, one can ob- 
tain the expression of the bare quasi-particle dispersion 
in terms of a bond charge quantity I^ — that 
in the momentum space is given by I^ = X]fc ^ka {~^) 
in a way that the kinetic part of the Hamiltonian 
©(neglecting the Stoner exchange) reads Hjj-mf — 
J2ko- [(1 ~ ^ji/ff - 2j2l-a)£k - rika- It is possible to 
show that the values of ji and j2 are linked to the V, J, 
and J' interactions by means of the following relations: 
ji = (J - V)/w and j2 = [J + J')/w[I3- Hence, by 
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FIG. 1: From the bottom to the top: superconducting gap, to- 
tal magnetization, and relative energy between the SF and the 
SC phase, respectively, as a function of a, for different values 
of the average effective mass renormalization. The density of 
electrons has been fixed to n = 1.2 and the couphng strength 
is g/t = 1. 

following the same procedure as in Model I, it is conve- 
nient to introduce the parameters a and 6, which now has 
to be determined self-consistently via the following rela- 
tions: 6 = [1 - (ji -I- j2)(/T + ^i)] > a = O2 - ji)(/T - ^i) 
Still, /cr depends on a and b by means of the occupation 
number of the a {(3) bands, via the following relations: 

Uk] = vl + ulinl + nl) - nl 
nui = vl+ul{nl+nl)-nl. 

where ul = (i)(l + -^=l^=pL=) and vl = I - ul are 

the coefficients of the Bogoliubov transformation. 

Few comments are worth mentioning on the possible 
solutions and on the effect of depaired electrons in the 
amplitude of the average mass b within the SF state. 
First of all, the amplitude of a and b is determined by 
j- = (j2 - ji) and j+ = {ji +32), respectively. More- 
over, the sign of j+ controls whether the starting value 
in SC state for b is larger or smaller than 1, while j_ sets 
the sign of the magnetization. As far as the stability is 
concerned, like in the Model I, it is crucial to see whether 
the introduction of depaired electrons in the SF phase 
reduces or increases the average effective mass with re- 
spect to that of the SC state. Hereafter, we will assume 
that j± are positive, which is consistent, in the weak 
coupling regime, with realistic values of the microscopic 
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couplings (J, J , V) above introduced. 

Let us now analyze the dependence of b in the case 
of the SF and SC phase. One can show that bsj gets 
undressed in the coexisting phase, thus its value is larger 
than bsc when 7^ 0. This result can be deduced by 
writing down the explicit expression of the ground state 
energy of the SF and SC phases. The expectation value 
of the bond charge term on the unpaired electrons gives 
a negative contribution to /| + /| which is not present 
in the non-magnetic SC case, thus increasing the average 
kinetic energy and consequently the value of Eq in the 
coexisting state. 

In Fig. [21it is reported the phase diagram for = 0.5, 
g/t = 0.25 obtained by varying with respect to the 
density. Of course, by modifying the amplitude of j+ one 
can span the phase diagram for all the possible values of 
ji and j2- We have checked that a change in j+ yields 
a shift in the critical line but does not give any qualita- 
tive change in the phase diagram. The phenomenology of 
the transition between the SC and SF state is the same 
of that observed for the Model I. Indeed, when one goes 
through the critical line, the magnetization has a jump 
to its possible maximal value which depends on A and 
/i, while the SC gap grows continuously from zero. For 
completeness, in FigI21it is also reported the line of tran- 
sition from a paramagnetic to a FM state in the case of 
absence of superconductivity. The shape of the region 
of stability in the — j+),7i] diagram can be under- 
stood by noticing that the kinetic exchange is given by 
kex — ^1 ~ = o, ^ and a — 2j_m(n — 1) (for the pure 
ferromagnetic case). Thus approaching the half- filling 
limit (n — > 1) one has that a — > and consequently 
the kinetic exchange goes to zero so that it is needed an 
infinite value of j_ to get a spin polarized state. This 
consideration explains the asymptotic behaviour of the 
critical line as we get close to the half-filling limit where 
it becomes more and more difficult to have a magnetized 
state because the difference in the bond charge occupa- 
tion for opposite spin becomes close to zero thus requiring 
high values of the coupling constants to create a charge 
unbalance. 

In conclusion, we have studied the occurrence of ferro- 
magnetism and s-wave singlet superconductivity within a 
single band model where the magnetic moments are due 
to a kinetic exchange mechanism, both for a generic and 
for a specific case with bond-charge coupling. It has been 
shown that the depaired electrons play a crucial role in 
the energy balance and that only when their dynamical 
effect is such to undress the effective mass of the carriers 
which participate in the pairing, then the SF phase can 
be stabilized. As far as the specific case of a bond-charge 
kinetic mechanism is concerned, we have shown that the 
phase diagram has a peculiar dependence on the density 
of carriers, and that only when the system is far from the 
limit of exact particle-hole symmetry, the SF phase can 
be obtained without going to very large couplings. 
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FIG. 2: Phase diagram relative to the Model II, showing the 
transition(fulI line) from the SF to the SC state, and that one 
from a paramagnetic to FM state when A = (dashed line), 
as the difference {j- — j+) and the density of electrons are 
varied. Due to the particle-hole symmetry, the part for n in 
the range [0, 1] is just symmetrically related. 
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